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Lorentz invariance is not necessarily an exact symme¬ 
try of nature. Though all the experimental results show 
that the laws of physics are Lorentz invariant, several ap¬ 
proaches to quantum gravity, for example, suggest that 
at very high energies Lorentz invariance might be bro¬ 
ken. Surprisingly, such violations, even if associated with 
Planck scale physics, may be constrained by present ob¬ 
servations. 

To discuss Lorentz invariance violation (LIV) we adopt 
the kinematic framework based on modified dispersion re¬ 
lations. That is, we alter the relation between the energy 
and the momentum of a given particle but keep intact the 
dynamics and the energy-momentum conservation law. 
Such modifications are indeed suggested by certain ap¬ 
proaches to Quantum Gravity Q. 

Modified dispersion relations manifest in different 
ways [2]. They may affect the propagation of photons 
by inducing dispersion and birefringence over long travel 
times. They could modify the threshold for known pro¬ 
cesses, as in 77 —> e + e~ and the GZK reaction §|. 
They might also induce new processes not allowed in 
QED such as photon decay (7 —> e + e _ ), photon split¬ 
ting (7 —> 77,777,...), and vacuum Cerenkov radiation 
(e —> e7). The non-observation of such anomalous ef¬ 
fects, then, constrains the assumed dispersion relations. 

We consider, following Ref. Q, the modified dispersion 
relation for the photon 


to the right and left circular polarizations of the pho¬ 
ton Q. It may be, however, that EFT does not describe 
the leading effects of Lorentz invariance violation. In 
any case, EFT is a dynamical assumption that goes be¬ 
yond the kinematic framework we are following. We will, 
therefore, neglect polarization dependence in the disper¬ 
sion relation and assume that all photons satisfy Eq. 0 
with £ positive. 

In standard QED the photon splitting process 7 —> 3j 
does not occur because, due to energy and momentum 
conservation, the three momenta must all be parallel and 
the amplitude vanishes in this configuration. However, 
such decay is allowed in the presence of the modified 
dispersion relation in Eq. 0 for ^ > 0. Since in this case 
the process 7 —> 3y has no threshold, its effectiveness in 
constraining £ depends solely on the rate. So far, only 
a crude guess of this rate exist in the literature. In this 
letter, we provide an estimate of the decay rate r(y —► 
3y) and show that the observation of high energy photons 
from astrophysical sources provides an upper bound on 
£ stronger than previously believed. 

To compute the rate of the decay 7 —> 3y we write 
the modified dispersion relation in terms of an “effective 
photon mass”, 


where E 7 and p 7 are the energy and momentum of the 
photon, the scale M = 10 19 GeV is an energy scale close 
to the Planck mass, and £ is a Lorentz invariance vio¬ 
lating parameter. This dispersion relation preserves ro¬ 
tation invariance but not boost invariance. Thus it can 
hold in only one reference frame, usually identified with 
the rest frame of the cosmic microwave background. 

Within an effective field theory (EFT) framework, the 
photon satisfies the polarization dependent dispersion re¬ 
lation E± = p 2 ± £(p 3 /M), where the subscripts ± refer 


so that the photon satisfies E 2 = p 2 + to 2 . 

Notice that all the scalar products of photon momenta 
in the process 7 —> 37, as well as all the scalar products 
of momenta and polarization vectors, vanish for to 7 = 0. 
The presence of a small m 7 ( ^ 0) changes the kinematics 
of the process, allowing the final photon momenta to have 
small components perpendicular to the initial momentum 
and proportional to mz/Ej. In addition, each photon ac¬ 
quires a longitudinal polarization vector ~ p^/m 7 . As 
a result, all the non-vanishing scalar products of photon 
momenta as well as all the scalar products of momenta 
and polarization vectors are of the order of the initial ef¬ 
fective photon mass. Thus, apart from the dependence 


2 


of the decay rate on 1 /£/, which remains explicitly in 
the final result, the important energy scale for the rate 
is given by the “effective mass” of the initial photon (the 
“effective photon masses” of the final photons are smaller 
and typically of the same order of magnitude of the ini¬ 
tial photons mass, due to their dependence on the photon 
momenta). 

We believe that this argument shows clearly that the 
relevant energy scale for photon decay is not the energy 
E 1 of the initial photon, but its “effective mass”. And 
due to the large suppression M -1 of this effective mass in 
Eq. ©, even energetic photons may have a small effective 
mass (say ro 7 ~ 10 _2 ^ 1 / 2 MeV for E ~ TeV). 

To proceed with our argument, let us define the “effec¬ 
tive initial photon gamma factor” 7 = m 7 /l? 7 , multiply 
and divide the standard equation defining the differen¬ 
tial decay rate by m 1 and extract a factor of inverse 7 
from it. The remaining equation is identical to the decay 
rate computed in the reference frame where the energy 
of the initial photon is m 7 (which is the “rest frame” 
of our “massive photon”), if we compute the amplitude 
from Lorentz invariant terms (there may be other non- 
Lorentz invariant terms but certainly the contribution of 
the Lorentz invariant terms should be there). In fact, 
if the squared amplitude |.M)j a | 2 term is invariant (here 
a stands for the initial state and (3 for the final state) 
we can change frames in the usual manner to compute 
the decay rate, since the other factors are 5 < ' 4 \p / _ 3 —p a ), 
still valid since we are assuming energy and momentum 
conservation, and the invariant volume of the final par¬ 
ticles, which would not be affected with respect to the 
standard case if we consider that only the initial photon 
has an “effective non-zero mass” (in any event the order 
of magnitude of the decay rate would not be changed by 
taking effective photon masses as defined above also for 
the final photons or only for the initial photon). 

We think that these arguments justify to estimate the 
decay rate of this “massive photon” in its rest frame and 
boost it -with the usual Lorentz factor- to the laboratory 
frame (which practically coincides with the rest frame of 
the CMB) to obtain the desired rate. This calculation 
accounts for the contributions to the decay rate of the 
Lorentz invariant terms in the Lagrangian (those present 
in the usual QED) although there may be other LIV 
terms in the Lagrangian |j], in particular those terms 
of dimension five which give origin to the modified dis¬ 
persion relations we consider. The contribution of such 
terms, however, seems to be smaller than the one we 
computed here J2j, fj] . 

The “effective photon mass” will need to be smaller 
than the electron mass (as shown below, to 2 < 4to 2 ) 
for the photon splitting to dominate the decay rate over 
the process 7 —> e + e^. In fact we will find that the 
“effective photon mass” is much smaller than the elec¬ 
tron mass. Thus the electron mass should not be ne¬ 
glected in the dimensional arguments given in Ref. [ 2 ). 


In this reference it is inferred that the decay rate must 
depend of a factor E 5 ^ 4 M~ 4 at photon energies E well 
above the electron mass. However, as we have shown, 
it is the effective mass of the photons which provide the 
relevant energy scale, and not the initial photon energy. 
We find here that the rate contains an additional fac¬ 
tor of (-E 7 /m e ) 8 £M -1 , which comes from the inverse 
gamma factor {m 1 /E 1 ) times a factor (m 7 /m e ) 8 due to 
the Euler-Heisenberg Lagrangian we use to describe the 
interactions among photons. This Lagrangian provides 
the effective 4-y coupling given by QED at energy scales 
smaller than the electron mass. Note that we are not 
modifying QED to actually include a photon mass in the 
Lagrangian. We use the modified kinematics induced by 
an “effective photon mass”, but do not modify the dy¬ 
namics. 

In order to gain confidence in the procedure we pro¬ 
pose, we will start by applying it to two processes for 
which rate calculations have been published, i.e. 7 —► 
e + e _ and p —> p"/ |j]. Let us first consider 7 —► e + e - . 
In the center of mass frame, the decay rate of a “mas¬ 
sive photon” into an electron-positron pair is simply 
Tqm — (a/2)m 7 . In the laboratory frame the rate gets 
an additional inverse gamma factor, 7 _1 = m 7 /I? 7 , 


_ a m 2 _ a £,E 2 

lab ~ 2 ~E^ ~ 2~M' 


( 3 ) 


This result coincides with that obtained in Ref. [fj, which 
is T ~ (a/2)(c 2 — 1 )E, for a photon with four momen¬ 
tum ( E,E/c ), if we replace the effective photon mass 
in Eq. (3) by its value in the model of Ref. f|, namely 
in 2 = E 2 (c 2 — 1). The energy threshold for this decay 
is obviously given by the condition m 2 = 4m 2 , which is 
equivalent to the condition E 2 = 4m 2 /(c 2 — 1) found in 
Ref. [I]. Thus, in terms of the “effective photon mass” 
the expression for the energy threshold and the decay 
rate are particularly simple. In fact, in the footnote 7 of 
Ref. 1,31], the procedure we follow is mentioned as provid¬ 
ing the right result. 

The rate for the p pj decay of Ref [ 6 | can also be 
obtained from kinematic considerations similar to those 
used above for 7 —> e + e _ . Let us use the modified dis¬ 
persion relation for the proton also proposed in Ref. 0, 


E 2 r = pl + m^ + 1 ^. (4) 

Since the effective mass of the proton depends on its mo¬ 
mentum, the mass of the parent proton, mp, is larger 
than that of the daughter p, nip, and the decay p —> py 
is kinematically allowed. In the rest frame of the parent 
proton the rate is Tqm — amp(l — rn 2 D /m 2 P ). In the lab¬ 
oratory frame the rate gets an additional inverse gamma 
factor, 7 _1 = mp/E p , thus 

T p ^ P7 ~ a— (1 - m D /m P ) ~ a— . (5) 
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The result quoted in Ref. 0 is for the rate of energy 
loss, dE/dx ~ —aEp(cp — 1) (see the notation of Ref |3( 
in which we have taken the maximal possible speed of a 
photon c 7 = 1 and the maximal speed of the proton to 
be c p > 1). Here c p = dE/dp , and from Ref. 0 c p ~ 
1 + r]E p /M. Since dE/dx ~ — TE p , the result of Ref. @ 
is the same we obtain. 

Let us now concentrate on 7 —> 3y. At energy scales 
(given by the “effective photon mass”) smaller than the 
electron mass, the interactions among photons are de¬ 
scribed by the Euler-Heisenberg Lagrangian, 


C-e-h = 


2a 2 
45 m 4 


_ p pnv 
2 ^ 


+ 7 ( -xcf lv pa F pv F pa 


( 6 ) 

Using this Lagrangian we estimate T(y —> 3y) in the 
laboratory frame (i.e. the CMB rest frame) as 


T(7->3y) = 


2 a 2 

~45~ 


1 


m® m 7 


3! 2 11 tt 9 ml E~ 


xf 


£ 5 E 

= 1.5 x 10 - 20 * 


M 5 m 3 


x / , 


( 7 ) 

( 8 ) 


where / is the integral over momenta with some large 
factors taken out, i.e. 


f=K 

m . 9 . 


d 3 ki d 3 k 2 d 3 k 3 4 


E\E2E 3 


5 A (p 1 -k 1 -k 2 -k 3 )\M\ , (9) 


and M is the invariant matrix element obtained from 
Eq.® by omitting the global factor 2 a 2 /(45m 4 ). 

Defined in this manner the factor / is of order one. 
Indeed, a similar (but not identical) integration was per¬ 
formed numerically in Ref. 0] and / turned out to be 
0.2. In our formulas we will keep explicitly the depen¬ 
dence on / and we will see that our final bound depends 
very weakly on / (in fact as / _1,/5 , thus only an / of order 
either 10 5 or 10 -5 would change the order of magnitude 
of our calculation). 

Thus the photon lifetime is 


t (7 —> 3q) = 0.025£ 5 f 1 sec. (10) 

For a given £ and a particular photon time of flight 
t , the condition r ~ t defines a critical value for the 
photon energy, E c . Due to the strong dependence of 
the lifetime on the energy, photons with energies above 
E c would quickly cascade down before reaching Earth 
from a distance ct and therefore would not be observed. 
Thus, photons reaching Earth from a distance ct should 
satisfy E 7 < E c or equivalently r > t. Since photons with 
E-y « 50TeV coming from the Crab nebula -10 13 seconds 
away- have been detected {§), we get the constraint 


£ < 1.2 x 10- 3 /- 1 / 5 


50TeV\ 2 
£7 ) 


( 11 ) 


This bound is considerably more restrictive than the 
bound £ < 10 4 previously obtained 12 ] using the same 
decay mode, and it is comparable to the strongest con¬ 
straint on £, i.e. |£| < 10 -4 9], coming from vacuum 
birefringence (different speeds for different photon polar¬ 
izations). This latter bound, however, was derived within 
the EFT framework, in which the Lorentz invariance vi¬ 
olating parameters for left and right circular polarized 
photons have opposite signs. Notice that the factor / 
that we do not compute in this paper would need to be 
~ lO -35 to bring the bound we obtain here to coincide 
with the bound mention in Ref. [ 3 ] , but since we have ex¬ 
tracted from the integral defining / all large dimensional 
quantities, we believe that the remaining integral cannot 
either very large or very small, but of order one. 

Notice that as soon as the decay 7 —> e + e _ is kinemat¬ 
ically allowed, i.e. for m 2 > 4m 2 , this decay mode, that 
happens at tree level, dominates over the one-loop sup¬ 
pressed 7 —> 3y We need to make sure that m 2 < 4m 2 , 
because then is the decay 7 —► 3y important. The con¬ 
dition m 2 < 4m 2 , which guarantees the validity of our 
approach, translates into a restrictive bound on £ less 
restrictive than the bound found above 


£ < 10 


-1 


50TeV 


( 12 ) 


Thus the required condition m 2 < 4m 2 is automatically 
fulfilled. This last condition insures also that the high 
energy photons coming from the Crab nebula did not 
decay into e + e“ before they reached Earth. 

In conclusion, we have estimated the decay rate of the 
process 7 —> 3y for photons fulfilling the Lorentz in¬ 
variance violating dispersion relation E 2 = p 2 + £p 3 /M 
and found that the observation of high energy photons 
from the Crab nebula sets an important constraint on the 
Lorentz invariance violation parameter £, much stronger 
than previously claimed. 
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